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ABSTRACT 


This thesis is a study of performance measures for 
acoustooptic spectrum analysis, and optical excision of 
narrowband interference. A Gaussian laser beam model is 
developed and programmed for propagation through thin lenses. 
Effects of beam truncation by Bragg cell apertures are ex- 
amined. Bragg cell performance is analyzed with respect to 
diffraction efficiency, information capacity, and response 
agility. A performance comparison is made between charge- 
coupled device and photodiode technologies in photodetector 
array applications. Overall excisor performance is discussed 
in terms of interference removal effectiveness and process 


corruption of residual information bearing signals. 
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1 ENTROBUETTON 


A. PROBLEM OVERVIEW 

The increased use of spread spectrum communication and 
ranging techniques has necessitated rethinking of intercep- 
tion scheme implementations. Two common requirements of many 
receiver systems operating in this wideband environment are: 
adequate spectrum analysis, and effective excision of narrow- 
band high power interference. Acoustooptic spectrum analysis 
and optical interference excision are two demonstrable methods 
of fulfilling these respective requirements. 

An acoustooptic spectrum analyzer transforms a wideband 
time domain signal into a spatially detected frequency domain 
signal. An Optical excisor at the detector can also utilize 
this space-time transformation. It performs real-time notch 
filtering of narrowband interference. Specific overational 
applications of the two devices include broadband receiver, 
radar and surveillance systems. All three military services 
rely on these systems and need good interference rejection. 
To meet this need, a NAVELEX 350/DARPA program was funded for 
the development of optical excision techniques. 

Thesteshnieaizebjjectives of this effort are to: 

l. Survey potential applications and analyze technical 

requirements imposed by each. 

Ee Ee optical excision configurations 


and components. 





3. Develop/procure major system components. 
4. Fabricate and measure alternative broadband excisor 
configurations to verify analytical results. 

5. Fabricate and field test a brassboard model. 

6. Design a preliminary advanced development model. 
The principal civilian contractor for the project is PROBE 
Systems, Inc. PROBE, assisted by the development efforts at 
the Naval Research Laboratory (NRL) and Naval Postgraduate 
School (NPS), was tasked with attainment of these objectives. 
NRL's specific program assistance was the development of a 
photodichroic optical clipper to be evaluated by PROBE. The 
research summarized in this thesis is part of the NPS con- 


tribution to the development program. 


PRO THESIS OBJECTIVE 

The objective of this work effort was to conduct a per- 
formance analysis of acoustooptic (a-0) spectrum analyzer 
and optical excisor components. The choice of this goal, 
over other possible goals related to the development program, 
allows maximal use of the Postgraduate School's substantial 
technical library and computer facilities. The approach was 
also felt to be the best way to augment PROBE's experimental 
study and minimize duplicate effort. This resulting study 
is intended to benefit the program by providing a comprehen- 
sive set of "first cut" performance guidelines. 

This objective was divided into four specific aims. Each 
effort relates to a particular system component. A summary 


of these goals is provided in Section E. However, more 
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complete explanations of a-o spectrum analysis and optical 


excision are required beforehand. 


C. ACOUSTOOPTIC SPECTRUM ANALYSIS 

An acoustooptic spectrum analyzer is primarily composed 
of an a-o beam deflector, photodetector array, and a laser. 
Its operating sequence begins with an input rf signal which 
excites a piezoelectric transducer mounted on the acoustic 
medium of the a-o deflector or Bragg cell. The transducer 
launches an acoustic wavetrain through the medium. The wave- 
train spatially modulates the refractive index and acts as a 
moving diffraction grating for the laser light. The acousto- 
optic interaction produces a one-dimensional angular pattern 
of diffracted beams. This pattern matches the spectral con- 
Ben or che original rf signal, Each individual beam in- 
@e@eity ts directly proportional to the corresponding rf 
signal frequency component intensity. Beam deflection angles 
are directly proportional to the component frequencies. A 
thin convex lens is used to focus the diffraction pattern on 
an image plane linear detector array. Lastly, the array 
photosites video-detect the beam pattern and transfer the 
feesieralWineormationtout electronically via shift registers. 

This method of wideband spectrum analysis offers many 
features which are difficult to achieve by other means. For 
example, a conventional scanning spectrum analyzer suffers 
from several shortcomings when used as a wideband surveillance 
receiver. Signals, which utilize direct sequence and frequency 


hopping spread spectrum modulation techniques, prove difficult 
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to adequately monitor. Pulsed signals are also troublesome 
because their probability of intercept is low. An a-o spec- 
trum analyzer readily accomplishes these surveillance tasks. 
Since the instantaneous bandwidth is very wide, probability 
of intercept problems are greatly diminished. The demonstra- 
ted output frequency resolution is excellent. These two 
features permit highly effective operation in densely popu- 
lated signal environments. Also, the Simpler physical com- 
ponents make the resulting equipment package smaller and less 
expensive than comparable conventional systems, especially in 
tne case of integrated optic system versions. 

Several articles which discuss currently available com- 
mercial units are listed in the bibliography. Reference 1 
describes an ITEK a-o spectrum analyzer with a 500 MHz band- 
End 1 MHz frequency resolution. Typical applications 
for this processor include broadband analysis and cueing 
receivers. 

Reference 2 discusses GTE Sylvania's work to develop 
acoustooptic wideband receiver and direction finder systems, 
presenting an angle of arrival detection scheme which utilizes 
multichannel Bragg cells. A GTE receiver having a l GHz band- 
width and frequency resolution of 2.6 MHz is also mentioned. 

A snapshot overview of acoustooptic development programs 
is related in Ref. 3. The article discusses recent advances 
made by the above two companies and related work done by 
others. It also examines exploratory work in integrated optic 
circuit chips. Hughes, Rockwell International, and Westing- 


house are the major companies deeply involved in this. 
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In addition to the above articles, the bibliography also 
notes a few landmark research papers. Reference 4 is the 
first article in this category. It is an excellent review 
of acoustooptic historical aspects which covers progress 
Pegenning with Brillouin’s 1922 prediction of light scattering 
by thermally generated acoustic waves. In this work, C. G. 
Quate and his associates neatly condensed much of the primarily 
academic work done prior to the 1960's. They also conducted 
a thorough two-dimensional analysis of wave interactions. 

The invention of the laser in the early sixties signifi- 
cantly increased the practical application possibilities of 
acoustooptic devices. Some of the first attempts to utilize 
the Bragg cell are discussed in Ref. 5. The authors, A. 
Korpel and others, explored the feasibility of acoustooptic 
television displays. They presented, in addition to the 
experimental findings, an analysis of one of the first a-o 
beam deflectors having a stepped transducer. The special 
transducer steered the acoustic wavetrain to maximize the 
proper light-sound interaction. This “beam steering” tech- 
nique increased the effective operating bandwidth of the 
Bragg Cell. 

R. Adler, one of the pioneers who worked with Korpel, 
gave a clear summary of beam deflector principles in Ref. 6. 
He also discussed a wide variety of acoustooptic applications 
which included television displays, spectrum analyzers, and 
acoustic imagers. 

Other references, which relate to specific analyzer 


aspects, are cited in following sections. However, the above 
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three papers best provide the flavor of development work and 


concisely present the significant acoustooptic principles. 


Pi OPTICAL EXCISION 

The optical excision process utilizes all a-o spectrum 
analyzer components with the addition of an optical or elec- 
tronic clipper. Only minor modifications are made in the 
operating sequence. The image plane diffraction pattern is 
generated in the same manner as previously discussed. The 
clipping action takes place in, or in front of, the image 
plane. Phase integrity is important for the time domain 
reconstruction of the signal, thus heterodyne-detection 
replaces video-detection to preserve signal phase. 

If an electronic clipper is used, then the array photo- 
Site outputs must be individually serviced. Outputs which 
correspond to frequencies in the notch band are ignored. The 
remaining outputs are summed to produce an excised version 
Sr che input signal. 

An optical clipper requires a single large area photo- 
detector to recover the heterodyned signal. An opaque ob- 
struction blocks or clips the diffracted beams produced by 
frequencies in the notch band. The remaining beams are 
mixed with a reference local oscillator beam and recovered 
by the detector. The net output is again an excised version 
EEN ie anal. A diagram of the optical clipping pro- 
cess is shown in Figure l. 

Optical excision has several advantages lacking in 


conventional approaches to narrowband interference rejection. 


14 





souenbas ssa0901g burddr12 testado `T oanbty 
indiano POSTOXH 


dNW 


/N 


SUETAd obeuy O 
Zeddrrz 
TeoTIdo y AOAATW 
SUIT] 
WIOFJSUPAL 


1192 Oberg 


AOAIATW 





r5 





AT n ci 
|—__——_—__—2l1mt2t17@€Ò@11ZzÀ1ÀZX2kpl1—_—————___—_—È—ÈÈ—È_n ; 


The process offers an arbitrary notch bandwidth with nearly 
infinite skirt slope. Optical excision lends itself well to 
adaptive notch filtering, which is important in signal environ- 
ments having time-varying interference sources of arbitrary 
number and center frequencies. 

Several efforts to develop suitable clippers are under- 
way. Photodichroic crystals are one set of possible con- 
tenders for application in optical clippers. These crystals, 
when used with crossed polarizers, act as optical attenuators. 
The amount of attenuation is dependent uvon the incident light 
intensity. This property gives them an inherent spatial fil- 
ter capability. Strong interference signals are sharply 
attenuated or clipped, and weaker signals of interest pass 
through. As mentioned in the overview, NRL is deeply invol- 
ved in this development area. Reference 7 provides a summary 
of experimental work on photodichroic spatial filters per- 

ol ed by WwW. ©. Collins,at the Naval Research Laboratory. 
PROM's or Pockels Readout Optical Modulators also show 
Ezawise in this area. A PROM optical clipper uses two lasers 

of different wavelengths. The first or "write" laser has a 
shorter wavelength and creates photoconductive spots on the 
PROM. These areas correspond to the interference frequencies 
to be excised. The PROM, positioned between crossed polar-- 
izers, is "read" by a laser of longer wavelength. Diffracted 
beams produced by the interference signal are blocked by the 
combined action of the photoconductive regions and crossed 
polarizers. The remaining diffracted beams are transmitted 


through the PROM and recovered by the heterodyned detector. 
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Reference 8 describes a PROM developed by ITEK which was used 
as part of an acoustooptic spectrum analyzer. 

The LCLV or Liquid Crystal Light Valve is another possible 
optical clipper. Hughes Research Laboratory has been the 
principal developer. The device consists of a liquid crystal 
photoconductor sandwich which uses scattering centers to dis- 
rupt the coherence of incident laser light. The scattering 
centers are created by the interaction of a cross-sandwich 
electric field and the incident light intensity pattern. The 
mesgente light lowers the resistivity of the photoconductor 
Samewach portion in proportion to thewintensity. This action 
increases the electric field in the liquid crystal layer and 
causes molecular alignment rotations. The areas of molecular 
rotation are the scattering centers. In order to be used as 
an optical clipper, the light valve is positioned between 
crossed polarizers. Attenuation is dependent upon molecular 
rotation and thereby the incident light intensity. Reference 
9 discusses progress made by Hughes on the LCLV. The article 
also compares various spatial filters' performance including 


the three mentioned here. 


Bee SPECIFIC GOALS 

The specific goals of this work were narrowed down to 
Four principal areas. First, laser beam profile effects on 
signal frequency component resolution required study. Second, 
Bragg cell performance factors needed determination for optimal 
system design. Third, photodetector capabilities and effects 


required scrutiny to compare the advantages of competing 
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detector technologies. And fourth, optical excision needed 
modeling to determine process effectiveness. The following 
paragraphs expand on these points. 

Laser light beam propagation was the first major aspect 
studied. This subject is important because the beam's cross- 
sectional behavior is a function of the propagation circum- 
stances. The beam profile significantly affects both a-o 
spectrum analyzer and optical excisor performance. 

For example, a wide beam is desired at the deflector 
aperture. This condition maximizes the acoustooptic inter- 
action. Conversely, a small light beam cross section is 
necessary at the image plane to discriminate the input 
signal's frequency components. 

Ray optics can not be used to determine beam profile 
behavior through the various system elements, as it does not 
adequately describe beam propagation characteristics. For- 
tunately, the complex variable equations, which do govern 
beam transmission through space and thin lenses, are avail- 
able; however, they are cumbersome to use. Simple methods 
of dealing with them were needed to increase their useful- 
ness as performance gauges. In particular, a hand-held cal- 
culator program was desirable to determine optimal lens 
placement and beam propagation indices. 

coco pagarion through the Bragg cell's finite 
dimension aperture affects cross-sectional behavior. The 
abrupt windowing action truncates the beam profile and 
causes sidelobes to appear. These sidelobes, if sufficiently 


large, contribute erroneous frequency components to the image 
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plane. Investigation of this phenomenon was required to 
determine the degree of performance degradation. 

The second thesis goal was to analyze factors affecting 
Bragg Cell performance, To pursue this goal, it was neces- 
sary to define good deflector performance. Uchida / 10 / 

did this concisely in terms of three criteria: high 
diffraction efficiency, large information capacity, and quick 
response. 

BEET action efficiency is the ratio of a diffracted light 
beam's intensity to that of the incident light beam. The 
input acoustic power is one of the components which determine 
G@eeeraction efficiency. A higher input acoustic power will 
give a greater diffraction efficiency, but only at a price. 
Increased acoustic power heats the interaction medium and 
creates troublesome temperature gradients. Resolution 
suffers and information capacity is lowered. Damping of 
internal reflections is made more difficult, causing an un- 
desirable persistence of signals which are no longer present. 
Deflectors which exhibit high diffraction efficiency with 
low acoustic power are very advantageous. 

Information capacity can be measured simply as the number 
of resolvable diffracted beams. It can also be influenced 
by geometric factors. Stretching the interaction medium's 
acoustic path length will increase the number of resolvable 
frequencies. However, this approach will also slow respon- 
siveness to changing signal environments. 

Simultaneous achievement of all three criteria, high 


HE su let ier etiiciene), large information capacity, and quick 
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response is not possible yet, but performance may be adequately 
tailored to a particular application by decreasing one un- 
necessary feature in order to increase a crucial one. 

Photodetector performance was another important subject 
for analysis. Detector performance criteria include: dynamic 
range, resolution, and information bandwidth. For spectrum 
analysis applications, there are two contending device types. 
Either a Charge-Coupled Device (CCD) detector array or a 
Charge-Coupled PhotoDiode (CCPD) detector array may be used. 
Both arrays are composed of a linear grouping of photosites. 
The photosites detect and integrate incident light intensities. 
Their response sensitivity is matched to the laser wavelength. 
Dynamic range is limited by photosite saturation illumination 
and dark current. Resolution is limited by the number and 
Wear Dacking density of the photosites. Information band- 
width is limited by electronic sampling considerations. 

ET op-icali excision applications which utilize an opti- 
cal clipper, a single large area detector is required. Po- 
tential device types include photodiodes and photomultiplier 
tubes. A pnotodiode may be adequate for a particular system 
if the information bandwidth of the signal is not large. 
Photomultiplier tubes are available for wideband requirements. 
Dynamic range limitations are similar to those for detector 
arrays. Resolution is determined by the optical clipper and 
the diffracted beam spot size. 

A number of references are available which highlight 
important aspects of the different detectors. CCD principles 


of operation are presented in Ref. ll. Reference 12 analyzes 
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photodiode array noise considerations in depth, Reference 
13 includes mathematical examinations of video and heterodyne- 
detection for photodiodes and photomultiplier tubes. 

Electronic and optical clipping methods comprised the 
final areas of analysis. Models of both techniques were 
developed. The first model addressed optical clipping. 

The blocking action of an optical clipper may be modeled 
to a first order approximation as a narrow obstacle in a 
coherent plane wave light field. The corresponding image 
plane diffraction pattern can be described in terms of 
Fresnel integrals. These integrals may be approximated by 
numerical methods. The resulting computer plots helo one to 
aen lize diffraction effects produced by an optical clipper. 

The next model is primarily associated with electronic 
clipping. It is intended to provide a measure of signal 
degradation which is induced by the excision process. This 
degradation which is caused by the removal of the desired 
signal's frequency components in the excisor's notch band. 
For Direct Sequence Spread Spectrum (DSSS) signals, this 
Besults im attenuation of the correlation peak. By compari- 
son, degradation effects for frequency hopped signals are not 
considered here as these signals do not derive any benefits 


on the excision process. 


meen RoPECTIVE 
The fundamental purpose of this introduction was to ex- 
plain the motivations for the work effort. Primary justifi- 


cation for the undertaking was the same as for the sponsoring 
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development program--a real obstacle to effective utilization 
of spread spectrum techniques must be confronted. The pro- 
blem of narrowband interference cannot reasonably be expected 
to diminish with only passage of time, more probably it will 
increase. 

Acoustoontics promises a solution to this Eero ence 
problem. The necessary technologies are maturing, and several 
available products ably handle missions similar to the in- 
tended ones discussed here. Therefore, this promised solu- 


ler is very likely to be fulfilled in the near future. 
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II, PERFORMANCE GUIDELINE ANALYSIS 


A. PARAMETERS OF NEGLIGIBLE EFFECT 

Several secondary effects are present in both the a-o 
Spectrum analyzer and optical excisor which were not con- 
Sidered as being germane to the thesis. Consequently, they 
Weee MOC examined in great detail. This section is a brief 
compilation of those effects. 

One performance affecting parameter outside the analysis 
realm of this thesis is input noise contributed by the re- 
ceiver system. No operation in either spectrum analysis or 
optical excision processors decreases white Gaussian noise. 
Above 50 MHz, atmospheric noise diminishes, and receiver noise 
figure becomes an important consideration. The intended appli- 
cations for both systems are assumed to be above 50 MHz. How- 
ever, the technigues dealing with receiver noise figure are 
well documented elsewhere and therefore are not addressed in 
tnis report. In all further analysis, the system noise in- 
duced prior to Bragg cell injection is assumed negligible. 

External vibration is another potential problem in both 
a-o spectrum analyzers and optical excisors. Component 
mounting vibrations can introduce undesirable signal modula- 
Mime Ait Currents may also disturb the signal. Coherent 
detection schemes are particularly sensitive to both occur- 
rences. Susceptibility to small ambient air currents was 


documented by M. King and his associates in Ref. 14. They 
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observed random phase modulation of a coherently detected 
signal which was caused by minor air movement. The experi- 
mental work of this thesis was done under controlled condi- 
tions with minimal mechanical disturbance. Also, video- 
detection, which reduced vibration concerns, was employed. 

In general, vibration considerations cannot be effectively 
encompassed within the framework of this report. The subject 
must be addressed with respect to each particular application. 

Laser frequency stabilization requirements of the optical 
excisor coherent detection scheme were also briefly investi- 
gated. Stabilization was concluded to be unnecessary for this 
application as the reference local oscillator is derived from 
the same laser beam as the information signal. Frequency 
Variations in tne local oscillator are matched by those in 
the signal beam. Only if there were a large path length dif- 
ference between the signal beam and local oscillator would 
there be any chance of beat frequencies appearing in the 
information band. This condition does not exist in the opti- 
cal excisors previously described. 

Tyndall scattering is a consideration in systems which 
employ liquid medium Bragg cells. It is not a concern with 
High purity crystal medium Bragg cells. Tyndall scattering 
is due to foreign suspended matter in the interaction medium. 
Scattering of this type contributes to background noise. It 
can be controlled by careful filtering of the liquid medium. 

Liquid medium Bragg cells are restricted to operation 
under 100 MHz due to attenuation considerations. Most new 
Bragg cells employ low acoustic velocity crystals. The a-o 
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EE EE used in experiments related to this thesis, housed 
a crystalline interaction medium, therefore, discussion of 
this topic is limited here. For those interested, further 
remarks on Tyndall scattering are made in Ref. 15. 

Temperature gradients in the Bragg cell medium have al- 
ready been described in the introduction as harmful to infor- 
mation capacity. Reference 10 documents experimental obser- 
maerons of thermal gradient effects. This condition can be 
controlled in essentially three ways. First, keeping the 
input acoustic power low minimizes input heat. Second, em- 
ploying materials with high thermal conductivity helps to 
Maintain uniform heat content throughout the medium. Third, 
id external cooling aids in dissipating heat introduced 
by the transducer. 

Demunıtorzmity of photosite responsivity for photodetector 
arrays is the last topic considered in this section. Typical 
variation for both CCD and CCPD arrays is ten percent above 
or below the mean. Normally, the intended applications do 
not require more exacting tolerances. In tne event that such 
precision is required, digital compensation techniques may be 
employed to offset individual photosite deviations. Reference 


ii Eontains more discussion on this subject. 


B. LASER BEAM PROPAGATION 

The first goal of the propagation study was to develop 
simple methods of dealing with the equations which govern 
transmission through space and thin lenses, and enable deter- 


mination of beam spot size and phase front radius anywhere 
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within the optic train path, In order to explain the results 
of the effort more concisely, a brief review of laser propa- 
gation fundamentals is presented, This review is based on a 
treatment provided in Ref. 16. 

The type of laser beam source discussed throughout this 
paper is a standard Helium-Neon unit which radiates at a 
.6328 micron wavelength. The beam has a Gaussian intensity 
profile and a spherical phase front. Beam spot size is de- 
fined as the perpendicular distance from the beam center- 
line to a point where the intensity has dropped by 1/e” Fictor: 
Beam radius is the phase front curvature radius. These two 
quantities uniquely determine the characteristics of free 
space Gaussian beam propagation. The plane perpendicular to 
the propagation path, where the spot size is a minimum and 
the radius is infinite, is referred to as the "waist" and is 
shown in Figure 2. 

Two equations determine beam spot size and radius as a 
function of directed propagation distance. Plots of beam 
radius and spot size, for propagation through free space, are 


depicted in Figs. 3 and 4. The equations are: 


Beam Radius = R(z) = z + (awe (1) 
o r f 2 2 ls 
Beam Spot Size = w(z) = w il + (1z/ (rw 1) } (2) 
where: 
ui” Minimum Beam Spot Size 
A = Laser Wavelength 


Propagation Distance from Waist 


N 
Il 
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Phase Front Radius in Decimeters 
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Figure 3. Phase Front Radius as a Function 
of Propagation Distance 
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Figure 4. Beam Spot Size as a Function 
of Propagation Distance 
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In order to predict beam behavior through a thin lens, 


the quantity q must be introduced. This value is known as 


ene” complex beam radius'. It is defined by the following 
relation. 
Complex Beam Radius = 1/g(z) = 1/R(z) - SAW” (2) (3) 


At the waist, the beam radius R(0) goes to infinity, leaving: 
do) = EEN = 4) 
g me 5 TE ( 


Reese manipulation of Equations (1) through (4) will yield 
Equation (5), which completely describes Gaussian beam trans- 


mission in free space. 
2) S= q, + 2 (5) 


Spherical waves which pass through thin lenses obey the 
lens law. For Gaussian beams, the complex beam radius under- 


goes the following transformation shown in Figure 5. 


1/q, = 1/4; - 1/f (6) 
where: 

q, = incident beam's complex radius 

q, = exit beam's complex radius 

f = thin lens focal length 


The above equations provided the basis for the TI-58 or 
TI-59 hand-calculator program to deal with Gaussian beam 
propagation through free space and thin lenses. This program 
is listed and explained in Appendix B. The first practical 
application of the program was to compute the beam spot size 


incident to the image plane photodetector array of an 











Figure 5. Gaussian Beam Transiting Thin Lens 
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experimental acoustooptic spectrum analyzer constructed as 
post of this work effort. 

The calculation's accuracy was dependent upon one pri- 
mary assumption. The beam profile could not be significantly 
truncated by the Bragg cell aperture or any of the thin lenses 
in the beam path. This assumption was adequately justified 
for the lenses used. A lens diameter of three times the spot 
ME is sufficient to contain 99 percent of the beam power. 
All lenses in the experimental set-up were greater than five 
spot sizes wide. 

However, this criterion was not upheld as the beam passed 
through the Bragg cell aperture. The beam spot size was 
expanded to completely illuminate the aperture and thereby 
enhance information capacity. Consequently, the beam pro- 
file was substantially truncated. In the limit, the incident 
beam cross-sectional intensity would have become uniform in 
distribution. The resulting exit beam, were it produced by 
the interaction of a plane sound wave and a rectangular light 
aperture of finite dimensions, would have had a horizontal 
intensity profile given below by Equation (7). This formula 


is based on the Fraunhofer approximation and is derived in 


Ber, 17. 
T(r,8) = I_(r) ¡sinfADO/A)) * (7) 
= PRO AAN 
where: 


Meee) = intensity aS a function of angle off the Bragg 
cell aperture boresight and path length to the 
image plane 
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Boresight light beam intensity 


D = Aperture horizontal dimension 
A = Laser wavelength 

6 = Angle off boresight 

E = Distance to transform lens 


The 1/8? beamwidth was determined from Equation(7) by an 


iterative algorithm. The solution is given by Equation (8). 
mie Function Le Beamwidth = SBW = 1.4rA/D (8) 


In order to determine the seriousness of violating this 
Criterion, beamwidth calculations were made for both extremes. 
The first assumed an incident Gaussian beam with no truncation. 
The second assumed an incident plane wave. Results of the 
calculations were compared with an experimentally determined 
value. 

The equipment used to make confirming observations of 
image plane beam spot size is block diagrammed in Figure 6. 
Component devices are described in the following paragraphs, 
and dimensional quantities needed by the propagation program 
are provided in Table I. 

The laser utilized for the measurements was a Spectra- 

FI sies 115. It produced a .5 milliWatt continuous wave 

output which had a TEM, spatial mode and random polarization. 
Power line ripple and random noise contributions were limited 
to one quarter percent of the output power. The longitudinal 


mode spacing was 550 MHz. 
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TABLE I. PROPAGATION PROGRAM INPUTS 


LASER 


EXT SEAM SPOT SIZE 


EXIT BEAM RADIUS 


WAVELENGTH 


SPATIAL FILTER 


DISTANCE FROM LASER 


LENS FOCAL LENGTH 


COLLIMATING LENS 


FOCAL LENGTH 


TRANSFORM LENS 


FOCAL LENGTH 


DISTANCE FROM COLLIMATING 
LENS 
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.44 mm 
INFINITY 


.6328 um 


14.8 mm 


STA 


.304 m 





A spatlal filter followed the laser in the operating 
train sequence. It removed intensity fluctuations which 
were due to light scattering from the beam profile. Dust 
on the laser windows and output mirrors is the major cause 
of this disturbance, although lens imperfection may be a con- 
tributing factor as well. 

The spatial filter performed its cleaning task by focus- 
ing the beam through a small pinhole. This focusing action 
caused the scattered light to spread out from the central 
beam. Light remaining in the central beam had a smoothed 
Gaussian profile. The pinhole passed the major portion of 
the central beam and removed the scattered light. Slight 
truncation of the beam by the pinhole produced an Airy Disk 
diffraction pattern. This phenomenon is discussed by Born 
and Wolf in Ref. 18. 

The outer rings of the Airy pattern were removed by pass- 
ing the central beam through a second aperture or iris. The 
net far field beam profile had a nearly ideal Gaussian shape. 
This beam was collimated by a thin lens. Doing so ensured 
complete illumination of the succeeding Bragg cell aperture. 

The particular Bragg cell employed in this set up was an 


Isomet (1205-1-2). It used Lead Molybdate (PbMoO as the 


4) 
interaction medium. The bandpass center frequency was 80 
MHz. This translated to a deflection angle of 7 milliradians. 
Nominal rf input power was one and a half Watts. 

An 80 MHz, fixed amplitude, continuous sine wave was chosen 


as the input drive signal. A pure sinusoidal tone was used so 


that no part of the measured beamwidth was due to the input 





signal. The deflected beam was assumed to be a function only 
of the interaction between the Gaussian light beam and the 
Bragg cell aperture's size and shape. Discussion of profile 
affecting light-sound interplay is delayed until the next 
section. 

When the Bragg cell was excited by the input signal, two 
light beams exited the device orifice. The first beam or 
"zero order diffraction beam" was not deflected. The only 
information carried by it was in the form of amplitude modu- 
lation. A stronger input signal would have meant a weaker 
undeflected beam. This beam was not useful and therefore 
terminated by a opaque zero order stop. The second beam was 
deflected, and represented first order diffraction. It was 
passed through the transform lens andsbreucht to a focus on 
the image plane linear detector array for width measurement. 

The transform lens played either one of two different 
roles which depended upon whether the diffracted beam had a 
Gaussian or Sinc function intensity profile. If the beam 
were truly Gaussian and still well collimated as it exited 
the Bragg cell, then the image plane spot size would have 
been insensitive to changes in the distance between Bragg 
cell and the transform lens. The important distance which 
needed precise calculation would have been that from the lens 
to the image plane. In the case of a well collimated Gaussian 
beam, this span is nearly equal to the lens focal length. 

Conversely, had the exit beam exhibited a Sinc function 
wiéé-i-sity profile, then the distance from the lens to the 


image plane would have been arbitrary. Such a beam would 
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have appeared to emanate from a point source at the Bragg 
cell aperture. A thin lens positioned a focal length away 
would have collimated the beam. The resulting spot size, 

as shown in Equation (7), would have been a function only 

of the distance to the transform lens, and not of the trans- 
form lens-image plane separation. 

The observed image plane spot size was sensitive to 
changes in either separation distance. This result implied 
the exit beam had retained characteristics of both profiles. 
The optimal transform lens position was located approximately 
one focal length away from the Bragg cell and a focal length 
from the image plane. This "best arrangement" is confirmed 
by Whitman and his associates in Ref. 19. 

The image plane sensor used to detect beamwidth was a 
(CCD121H) linear detector array developed by Fairchild. It 
consists of 1728 sensor elements and is a charge coupled, 
buried channel device. The output is delivered serially 
and, for the experiment, went to a Tektronics 422 oscillo- 
scope display. A typical output display, which resulted from 
an 80 MHZ monotone input signal, is shown in Figure 7. 

The physical apportionment of the photo elements is 
Shown in Figure 3, The ie beamwidth was measured by ex- 
panding the oscilloscope time display and counting photosite 
output peaks within the beamwidth. This number was converted 
into a spatial dimension by using the known dimensions in 


Eidure 8, 
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The Sinc profile was determined using Equation (8) and 


the following values: 


r = focal length of transform lens = .304 m 
À = Hellum-Neon Laser wavelength = .6328 um 
D = effective aperture width = ./ cm 


The Gaussian beamwidth was calculated by using the propaga- 
[men program and Table I. Intermediate results are listed 
Mi Table 11. Incidently, this table may be used as a 
standard solution set to ensure all proper entries are keyed 
into the calculator's program memory. 

Table III compares the three independently determined 
values. These results indicate the actual beam profile was 
very Similar to the Sinc function. However, the highest 
sidelobe was measured 17 dB down from the central peak. If 
the profile were an ideal Sinc function, then the highest 
sidelobe would have been approximately 14 dB down. 

D. L. Hecht presented a concise sidelobe analysis of 
truncated Gaussian beams in Ref. 20. In one graph he related 
the maximum sidelobe level to beam truncation ratio- D/2w. 
This graph is shown in Figure 9. Again, this does not take 
into account effects of the sound column of the light beam 
ME Eile, it is a function only of truneation ratio. 

Figure 9 provideä a rough means of indirectly determining 
truncation ratio. For a -17 dB sidelobe, the corresponding 
eeben ratio was in the range of .55 to .6. The actual 
truncation ratio value was calculated as .47. It was deter- 


mined by using the spot size given in Table II for the beam 
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TABLE II. BEAM PROFILE INTERMEDIATE RESULTS 


At Laser Aperture 


Beam Spot Size w= .44 mm 


Beam Radius = R 


Il 
8 


After propagating .l m 


„4424 mm 
IT m 


W 
R 


At Spatial Filter Lens Exit 


.4424 mm 
-14.82 mm 


wW 
R 


After propagating 14.820044 mm 


6.749 um 
R = 0 


W 


After propagating .254 m 
w = 7.581 mm 
R = .2540 


eec ollimating Lens Exit 
ASSL mm 
ER mn 


WwW 


After propagating .3 m 
w = 7.581 mm 
R = 1.770 x 10° m 


Fi DE arstorm Lens Exit 
7.581 mm 
-.3040 m 


WwW 


Beem propagating 3.039997 m 
w = 8.077 um 
R = © 
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TABLE III. COMPARISON OF BEAMVIDTH VALUES 


Sinc Image Plane Beamwidth 38 um 
Measured Image Plane Beamwidth 34 um 
Gaussian Image Plane Beamwidth l6 um 
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entering the Bragg cell (.75 cm) and the deflector aperture 
width (.7 cm). The graph determined value differed by 17 
percent. 

This disagreement is not very surprising in view of pos- 
sible sound column effects and other error sources. Channel 
stops were believed to have played a role in limiting pre- 
cision measurements. The intervening metal strip between 
two adjacent photosites consumes 24 percent of the available 
surface area. Its effect can be minimized but not eliminated 
Ey rocking” the diffracted beam back and forth. 

Maen On a precision focus could have also distorted ob- 
servations. If present, minor lens and Bragg cell imperfec- 
tions would have slightly affected ultimate focusability of 
the beam, however, operator error was not particularly sus- 
pect here. Even though focusing was a tedious exercise, 
minor lens vosition changes produced very noticeable amplitude 
deviations in an expanded scale oscilloscope display. 

Also, the detector array was not a perfectly one-dimen- 
Sional sensor. Each pixel averaged the incident light in- 
tensity over a rectangular area. The vertical photosite 
dimension was twice the magnitude of the horizontal. It was 
also equal to the main lobe beamwidth calculated for an ideal 
Gaussian beam. It is not fair to assume that the ratio of 
averaged intensities would have necessarily been identical to 
the ratio of the peak values. A closer ratio agreement would 
have been realized with a smaller vertical dimension. 

ile subject of diffraction beam sidelobe effects is the 


last remaining topic of this section. D. L. Hecht / 22 / has 
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shown a way to anticipate sidelobe levels, but no one, to this 
author's knowledge, has specifically addressed sidelobe sup- 
pression. The obvious approach of increasing the truncation 
ratio to suppress sidelobes has the previously mentioned 
drawback of reduced information capacity. In present systems 
having only 30 dB overall dynamic range, a -20 dB maximum side- 
lobe level may not be significantly harmful. However, wider 
dynamic ranges in newer systems are now evident, so eventually 
the problem must be confronted. 

Since spatial sidelobes are a direct result of the abrupt 
windowing action by the aperture, the only way to remove them 
is by making the absorption-transmission changeover a less 
discontinuous process. To produce a smoother window function, 
a graded film mask could be deposited on the deflector aper- 
ture exit. Or, absorptive impurities could be diffused into 
the interaction medium to produce the same effect. The ab- 
Fi pion profile would only need to be a function of the hori- 
zontal axis as vertical sidelobes do not interfere normally in 
a one dimensional array. Figure 10 depicts a possible graded 
transmission profile which exhibits a smooth transition from 
Bnzinite to zero attenuation. 

WpisW$compilletes the formal discussion on laser beam pro- 
pagation. To recapitulate, a brief review of propagation 
characteristics was presented. Behavior-determining equations 
explained in the review were utilized in a calculator program 
for beam spot size and radius prediction. Gaussian beam model 
results were compared with experimentally observed figures and 


those predicted by a plane wave-aperture interaction model. 
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Truncation ratio was found to play a dominant role. Finally, 
the significance of sidelobe interference was discussed and 


two methods of sidelobe suppression were suggested. 


C. BRAGG CELL PERFORMANCE 

As mentioned in the introduction, a Bragg cell translates 
a time domain signal into the frequency domain by converting 
a Signal's component frequencies into beam deflection angles, 
and respective intensities into beam intensities. This process 
is sketched in Figure 11. How well the Bragg cell accomplishes 
this task is the subject of this section. The three major 
performance criteria of diffraction efficiency, information 
capacity and response time were used to partition the investi- 
gation effort. The first topic under analysis was diffraction 
Siti Clency. 

Before beginning this discussion, a brief mention of the 
basic law governing Bragg deflection is required. This is 
given by Equation (9) which relates the deflection angle to 


the acoustic wavelength of the input signal. 


À 
geg, EB (9) 
21 
O 
or 
MN Af 
"29 
where: 
do = laser wavelength 
N = acoustic wavelength 
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6 = deflection angle 
f = input frequency 
V = acoustic velocity 


A detailed geometric proof of Equation (9) is given by Yariv 
Mi Ref. 13. Since 8 is very small, typically on the order 

of seven milliradians, the first order approximation is 
sufficiently accurate for spectrum analyzer and optical ex- 
cisor applications. Only the two factors of laser wavelength 
and acoustic velocity are able to alter the ratio of deflec- 
tion angle to input frequency. The laser wavelength is chosen 
from those available to best match the transmittance of the 
raso cell, and the acoustic velocity is a function of other 
constraints, as will be shown later. 

On the other hand, the beam diffraction efficiency has 
many factors which affect its value. This parameter is de- 
fined as the ratio of the diffracted beam intensity to the 
intensity of the original undeflected beam. 

This is shown in Equation (10) which determines the ratio. 


Mme was derived by I. C. Chang in Ref. 21. 


I] BEE {n + Ree 
een (10) 
O in + (AK,L/2) } 
with 
2 6 2 PL 
aa LEO y ‘Glas 
2 9 
25 oV Hcos 8, 


50 





u - diffracted beam intensity 
= Original undeflected beam intensity 
Ao = free space laser wavelength 


= refractive index 
p = elastooptic coefficient 


= mass density 


V = acoustic wave velocity 
2, = acoustic power 
L = interaction length 
H = acoustic beam height 
AK, = wavevector mismatch 
SR = deflection angle determined from equation (9) 


An approximate version easier to work with is given by 


Equation (12): 


Ne 
I EE) 
Een (12) 
O {AK,L/2} 
Let 

AS BEES {n3} 

ae” 

Sa {AK,L/2} 

B = EER EE EE 

(AK, L/ 2)” 


Term A of eguation (12) is constant with respect to input 
signal frequency. Its effect will be examined later. Term B 


fixes the relative bandshape of a normal acoustooptic beam 








deflector. By converting wavevectors to frequencies, B may 
be expressed as: 
sin°{n/2L/L_{FF_-P°}} 


(m/2L/L_(FF, -F N 


where: 
F = normalized center frequency = 1 
I = input frequency normalized to the center frequency 
value 
2 SE D 
L_ = nA cos /A_ AraC teristic interaction length 


ts teranstormation is shown by D. L. Hecht in Ref. 20. 
Deckert Of this equation with various ratios of L/L. is shown 
Mmeerigure 12. The 3 dB bandwidth is given by equation (14) 


and is also from / 20 7. 


Af = 1.8f _L_/L (14) 
where: 
N = center frequency 


Equations (10) through (13) describe Bragg regime or 
"thick grating" diffraction. Only a first order beam will be 
deldl'ected with a thick diffraction grating. A thin diffrac- 
tion grating will yield many higher order diffraction beams 
ii addition. This type of diffraction is referred to as 
Raman-Nath diffraction. It is intolerable for the intended 
applications because of interference caused by the additional 
beams. 

Equation (14) implies that the deflector bandwidth can be 


easily increased by shortening the interaction length. 
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Figure 12. Diffraction Efficiency Bandshape 
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Increasing bandwidth in this manner is of limited overall 
benefit as the diffraction shifts from the Bragg regime to 
Raman-Nath. In general, normal Bragg diffraction is restric- 
ted to one octave. 

This leaves only one effective way in Equation (14) to 
emearge the Bragg cell's operating bandwidth. The input 
signal's center frequency, Los must be increased. The ability 
to increase f is determined by the input signal transducer. 
Substantial transducer analysis is included in Refs. 10 and 21. 
The highlights are briefly presented here. 

The acoustic transducer consists of a thin piezoelectric 
strip bonded to the interaction medium. Its effectiveness in 
launching an acoustic wave into the interaction medium is a 
function of coupling and impedence matching characteristics. 
The particular bonding technique bounds the maximum frequency 
that the transducer can effectively couple to the medium. 

The degree of electrical and mechanical impedence matching 
Metermines the transfer efficiency from electrical to acousti- 
Bae signal energy. It is also frequency dependent. Reference 
2 reports some available deflectors which operate in the 

three to four Gigahertz region. 

The above discussion of tranducer effects completes the 
analysis of Equation (12)'s term B. Attention is now directed 
to term A. It relates the interaction medium's material con- 


memes to diffraction efficiency. 
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The argument of term A may be rewritten as: 


pe TT. 
n => C (15) 
where: 
c= 12 E E (16) 
N Zë 2 H così, 


Since n° is a sine function argument, improvement in diffrac- 
tion efficiency is a "simple" matter of increasing C. It is 
desirable to increase the refractive index - n and elastooptic 
coefficient - p, aná decrease the acoustic beam height H, 
the acoustic velocity - v and the mass density - p. Little 
can be done with Mo Y2 or 90° plico he deleterious efrects 
of increasing acoustic input power have been noted in the 
BEE Eiere And, inordinate increasing of L would prove 
costly in terms of bandwidth by Equation (14). 

Several acoustooptic figures of merit have gained accep- 


tance for use in the selection of suitable interaction mediums. 





Maey are: 
Di 2 2 
= — 7 
M, Se M, {nV } (17) 
> 2 
M, = —- (18) 
oV 
Al (19) 
RS = M, {nV} 
oV 
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M, = EI = M, {nv}? (20) 


Gordon introduced H. and M, in Ref, 22. His results were 
Et ther clarified by Chang in Ref. 21. Chang also included 
a discussion of M, and M,- My in equation (21) related diffrac- 


tion efficiency to acoustic power for a given bandwidth. This 


Bemyielded by Equations (11), (14) and (17). 





M 
n ~ 1 9 
P AE er v (21) 
O O 
where: 
9 = 1.812 


This same efficiency-power ratio can be expressed in terms 


of a fixed geometry by using M, in Equation (21). 


2 
n = L f 5 
D_ M {7 } i (22) 
a to 


M, is tne preferred figure of merit for systems in which the 
transducer height is not limited by impedence or fabrication 
constraints. H would then be limited only by the laser beam 


height entering the cell. For convenience, H is replaced by 


meet, Which transforms Equation (21) into Equation (23). 





a} (23) 
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where: 
Pa = P_/LH (25) 


In addition to having high values of the above figures of 
merit, a good interaction medium should exhibit other charac- 
teristics as well. Low coefficients of optical and acoustic 
attenuation are desirable. Optical attenuation degrades over- 
all performance and ocntributes heat to the material. Acoustic 
attenuation decreases information capacity. Lastly, if the 
material is crystalline, then consideration must be given to 
how easily a continuous crystal structure can be pulled from a 
melt. For specific information on medium performance, see 
Ref. 10 which contains an extensive comparison of material 
properties. 

Once the deflector material has been chosen and the device 
dimensions fixed, nothing remains free in Equation (10) to 
improve diffraction efficiency except for the acoustic input 
power. However, beam steering and other compensation techni- 
lies can still be called upon to improve performance more. 
Korpel's work / 5 / on stepped transducer beam steering has 
already been briefly cited. Gordon / 22 / also discussed beam 
steering, only with a planar instead of stepped array. Both 
transducers represented first order steering as each successive 
element differed in phase by m radians. 

Pinnow Ni compared normal or single transducer element 
Bragg diffraction with stepped and planar first order beam 
steering. He argued effectively that beam steering was better 


as it required less total acoustic power and a lower cell 
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acoustic power density. This is shown most clearly by attemp- 
ting to double bandwidth in a single element deflector cell 
by doubling center frequency. The effects are revealed in 


Equation (14) with L_ rewritten in terms of frequency. 


2 2 

E 2 _ n V ” V 

Lo = nA¿cos8 A3. = i Bee = BS (26) 
O 


“2 
——) (27) 
e, 


f = 1.8£ L/L = 1.8f KI 
Boubling the bandwidth by doubling Li causes the transducer 


length to reduce to one quarter as show in Equation (28). 


2 
V 
e! (28) 


2Af = {1.8}{2f K}{ 
2 4£°L/4 


This occurrence induces the acoustic power from Equation (22) 


Ho quadruple. 


m ta SH 
Ae 


O 


D 


The net result is four times the original acoustic power is 
Spread over one fourth of the original surface area. Pinnow 
explained how this unhappy conclusion could be largely cir- 
cumvented by using a first order array of only moderate com- 
plexity. 

e a me pointed out that a stepped transducer array 
is nearly twice as efficient as a planar array. A first order 
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planar array must divide its power between two beams which 

are symmetric about the normal. The corresponding stepped 
array utilizes its inherent asymmetry to put a large percen- 
tage of the acoustic power into a single main lobe. A typical 
first order stepped array is depicted in Figure 13. 

The performance improvement provided by a stepped array 
was Calculated in terms of an effective acoustic power ratio. 
To accomplish this, the array beam pattern needed to be 
Metermined. For simplicity, it was described in terms of far 
field approximations even though the light-sound interaction 
actually takes place in the near field. The overall pattern 
was taken to be the produce of the single element and array 


functions. 


ME) we) = B(8)C(8) (30) 


where B(9) is defined as the single element function and is 


given by 
ein E 2 
B($) = Sor: (3L) 
Ii } 


Mie) is the array function and is 


SS 2 (8-91) Se 
C(8) =. ——- 5 
NS cr 2 
A az 
9' = beam steering angle 
6 = angle from normal to element plane 
D = center to center element spacing 
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9 = o or ı phase shift 


Eie 1 First Order Stepped Array 
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Il 


element length 
N = number of transducer elements 

The array performance was compared at midband with a 
deflector having one element of length Lo by taking the ratio 
a of the effective array beam power over the effective single 


element transducer beam power. 


1 
P. LW 
O a (33) 
P LW 
a OS 


where: 
P, = Array Acoustic Power per unit lingth 
' 
= P_/N 
P = Unit Transducer Beam Power per unit length 
ne = Array Pattern value at midband 
We = Single Element Pattern at midband = l 


This leaves 





wot Wa Ha 
ola es 3 
N Lo N? 
where: 
g = NL/L_ 


is tne ratio of the array length to tne single element trans- 
ducer lengths. 
ee midband: 


Be = G< 1 (35) 
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> EECHER S 
N a (36) 
sintq-10-0')) 


C(8) 


= N? 
therefore: 
a = GN 037) 


The end result is a performance factor a equal to something 
less than the number of transducer elements. Pinnow did caution 
however, that increasing N much beyond four was self-defeating 
as the inherent errors of first order steering inhibit further 
improvement. Thus, first order beam steering can be utilized 
to provide some bandwidth improvement without increasing power. 
Its benefits must be subjectively weighed against the diffi- 
culties caused by increased device complexity. 

Although beam steering does somewhat improve diffraction 
ratio consistency across the passband, there are also other 
optical and electrical techniques which can further minimize 
nonuniformities. A set of staggered bandpass filters inserted 
prior to the deflector could pre-distort the input signal to 
even out the humped midband response of a normal Bragg cell. 
This method would be amenable to both spectrum analysis and 
optical excision systems. If advantageous, an optical excisor 
could also employ post-distortion filtering equally as well, 
however, a spectrum analyzer would required a scheme of time 
Varied gain for post-distortion. An optical filter or atten- 


uator whose transmission is inversely related to the diffraction 





ratio bandshape would be another possible contender. One 
method of obtaining such a filter employs photographic film 
as the filter medium. The technique tailors the component to 
a particular system by sweeping a single rf signal across the 
deflector passband while the resulting diffraction beam ex- 
poses an undeveloped film slide positioned between the de- 
flector cell and photodetector. This film, once developed, 
would be a post-distortion filter and provide a uniform pass- 
band response, as portions of the film which were exposed 
to higher illumination intensities prior to development would 
attenuate the diffraction beam to a greater extent. 

The analysis up to this point has assumed the deflector 
treats each rf signal frequency component independently of 
all others. This supposition is not entirely true as many 
secondary effects may be noticeably present. Each element 
of the composite input signal depletes the incident light 
beam's capacity to respond linearly to the other simultaneously 
oresent components. Diffracted beams may also be rediffracted 
by other acoustic wavefronts encountered prior to exiting the 
cell. D. L. Hecht / 24 / analyzed multiple frequency component 
effects. In order to do this, he used coupled mode equations 
to solve the classical wave equation, Equation (38), for N 
Signals in both the Raman-Nath and Bragg regimes. 


oA tal 3 E (38) 


fe 37 
C m 


His solution algorithm is similar to that followed by Klein and 


Cook in Ref, 25. They, however, only treated the interaction 


A 





of a single frequency sound column with the incident light 
beam. ee. development and results for Bragg diffraction 
are discussed here as this type of interaction is more impor- 
tant for the intended applications. 

The electric field E was represented as an N-dimensional 
Fourier series expansion which was much simplified for Bragg 
regime diffraction since only two principle modes were present. 
This allowed the single second order differential equation 


expressed in Equation (38) to be broken up into two first 


order coupled mode differential equations. 


Ay Oz a Von = 
qz Y N) = z 57 Y (nta) (39) 
m=1 
N V 
SÉ m u mi ed E 
a (n) = E 35 Y (n a) (40) 
mel 
where: 


Be aulescent index of refraction 


yO = zero order diffraction mode 
oe sR is. order diffraction mode 
Ze (My Mar... ny) = dimensions of the N-tuvle Fourier 


expansion of E where Dj Da: etc. 
have integer values between -~ and 


coe 
(n-a_)= GEET EE My By? 


Mia) (nnar e Dan Pmt] Pmt Ar Dal 
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ees at Conditions at z = 0 


I 


v° (0) 


tn) 


0 


Hecht's results to these equations are given in terms of Va 
which is "the normalized index of refraction amplitude" for 


the mth signal. 


where: 
e = Ku L/cos8 | (41) 
ue = refractive index perturbation due to mechanical strain 
caused by the mth signal. 
K = T/A, 


vm MA 1 ited to the diffraction ratio efficiency quantity 


by the following brief development. 


3 
rc EOS 42 
o A El n ERC (42) 
where: 
21, 
SC = = mechanical strain due to mth signal 
oV 


(Note that V is the sound velocity) 


In = Acoustic Intensity due to mth signal 
= P /LH 
Then 
6 2 

V u P L 

M 2 ME (43) 
2 8H cos 9 
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Equations (ll) and (43) allow us to equate optical diffraction 
efficiency We with a, Differentiating and decoupling equations 
(39) and (40) enabled them to be solved with assumed power 
series solutions. The results for N signals were given in 
terms of truncated series involving sé 

The diffraction efficiency of a single frequency component 
signal in the presence of N-1 other single frequency component 


signals was found to be the following when evaluated at Z = L. 


E V Wien 2 N Ee 
eae eee 1 N 7 ms i 
a (y (a,)} [>=] = cia 3 E Ge } (44) 
O 1=2 
D 
= {>} provided all V, < .l 


The corresponding fractional deviation from ideal diffraction 
ratio of one signal due to N signals has been defined as 


kompression!’ C. 


V 2 
o {=} +2 2. ES (45) 
1 2 3 2 A 2 
N 12 
13 
2 


e e à e ` 
'Cross modulation', m, is the compression of one signal due to 


a second signal only. 
M, 5 = $i} (46) 


¡Sede compression is not included here.) 





Compression and cross modulation are important considerations 
in high density signal environments. Cross modulation is par- 
ticularly so for optical excision applications, as an amplitude 
modulated interference signal of high power could do much harm 
to the reconstructed signal of interest. 

Intermodulation is another area of concern which is described 
ma terms of a coupled mode solution yo (a,-a,). Second order 
intermodulation products are the highest power intermods and 
contribute to the zero order diffraction mode. As long as 
, then the second order intermod products will not 


nin 


Be of concern. This is guaranteed if 2f . >f ; 
min max 


where: 


ER: 
min 


minimum deflector passband frequency 

Zu. = maximum deflector vassband frequency 
This condition speaks well for deflectors having a bandwidth 
of slightly less than an octave. However, two-tone third- 
Ber |itermods will sti fall within the deflector pass- 


band and produce first order diffraction beams whose duras. Lon 


efficiency is described in terms of p? (2a}-a,) as follows: 
"AE EEN EE (47) 


Three-tone third-order intermods of the form (£,+£,-f3) will 
also create first order diffraction beams in accordance with: 
y 2 y 2 2 


V 
ì Ly di (48) 


Dee’ 


Ol Fe 


(a,ta,-a3) Š 
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As system dynamic ranges improve, third order intermodulation 
products will continue to become more significant. 

These last remarks complete discussion on multiple fre- 
quency performance effects. Attention is now drawn to another 
closely related subject concerning geometric implications of 
light and sound interaction. Uchida / 10 /, in his derivation 
of diffraction efficiency, assumed weak interaction of plane 
wave light incident to a rectangular sound column. The re- 
sults, while very useful, do not take into account true geo- 
metric considerations of the light and sound beams. Several 
authors have challenged various aspects of the problem. Two 
particularly illuminating papers are briefly discussed here. 
The first paper / 26 / presented closed-form or analytic solu- 
tims to problems where the light and sound beams were of 
Similar profiles. 

The most common type of light-sound interaction experl- 
mentally observed is between a Gaussian profile light beam 
and a rectangular sound column. McMahon / 26 / argued that 
a closed-form description of this instance could not De 
arrived at. Instead, he proposed to bracket and indirectly 
analyze this case by finding solutions for two closely related 
cases. The two he chose were a light beam of rectangular 
profile intersecting a like sound column and a Gaussian to 
Gaussian intersection. He conducted the analysis by using 
an oscillating electric dipole model formulation and was able 
to determine the first order beam power for both situations. 
The following results were derived for the light beam incident 


at the Bragg angle so that wavevector mismatch would not have 


6 8 





Ro be considered. 


For the rectangular to rectangular interaction 


= a 
ee = Atl 3) (49) 
where 
SS = power in first order diffraction beam for rectan- 


gular to rectangular interaction 


2 
WL 


a G £(P,) eos.G 
O 


£(P_)= function of input acoustic power and material 


constants 
W = laser beamwidth 
E = sound column width 
_ 2Lsin® 


W 


The Gaussian-Gaussian interaction gives 


P =n *A/{1 ee (50) 
gg 
nere: 
E, = power in first order diffraction beam for Gaussian 


to Gaussian interaction 
Wiefeutnorclalimed for the condition of (*<a<2) that the 
power diffracted by either configuration would vary from the 
other by no more than eight percent. Thus, by taking the mean 
value between them, one could only be off by about four percent 


from the "true" value for the intermediate Gaussian to 








rectangular interaction case. However, for the deflector used 
in the experimental set up of this thesis, the calculated 
figures for the two configurations differed by a considerable 


margin since o was not within the specified limits. 


mo 


wr~ LZ 

de = 7 milliradians 
a = .014 
D = (.995)A 

rr 

P = (9.86)A 

gg 


In spite of this variation, the results do give some insight 
Bemgeometric influences on diffraction efficiency. To obtain 
a more general appraisal of the interaction process, Korpel's 
work in Ref. 27 should be consulted. He presents a technique 
for handling light and sound beams of arbitrary profile. The 
ensuing derivation assumed a spatially varying plane wave 
angular spectrum representation of the radiated fields and a 
"multiple scattering model" for the interaction process. 
Numerical results for particular beam configurations require 
power series solutions that the author left for future deter- 
Mination. 

In conclusion, primarily qualitative benefits are available 
from McMahon's analysis of geometric factors. And while Korpel's 
theory offers exact solutions for general beam profiles, im- 
plementation may be difficult in specific cases. Much further 
work remains in this area. 

The next and final sub-topic related to diffraction effi- 


ciency is a discussion of Bragg cell dynamic range. This 
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performance measure is defined as the linear operating range 
limited by optical background noise for minimum detectable 
signals and saturation for large signals. There are three 
types of optical scattering which are causes for background 
noise. Airborne dust and lens defects are scattering centers 
Wwielbwsdedtin the first class. Their noise contribution is 
basically a function of the environment and component quality 
for a particular operating system, and so will not be consi- 
Ed here. By comparison, optical scattering resulting from 
the diffraction of the zero order beam is a function of design 
and is one of the two possible theoretical limits for minimum 
background noise. The other possible limit is scattering 

due to thermally generated acoustic waves in the Bragg cell 
medium. This type of scattering in known as Brillouin 
scattering. The saturation limit for large signals is deter- 
mined by zero order depletion considerations in multiple 
frequency environments. Another maximum linear response limit 
could be used in signal environments composed of only one 
frequency component at a time, however this limit would have 
no practical value for typical operating conditions. 

To determine a system's minimum detectable signal, both 
theoretical limits must be calculated and compared for the 
larger value. The minimum detectable signal due to zero order 
diffraction effects is considered first. In the worst case 
(maximum sidelobes), the diffraction pattern of the main beam 


obeys Equation (7). Sidelobe maxima are roughly determined 


EE: 


TI 








T (0) TX (31) 
max x= (N+) 
- 2 
2 g 
TN = 


where: 


o 
Il 


D8/A, 


A 
Il 


angle of deflection due to the minimum passband fre- 
quency/angular width of a single spot and is given 
Dy ETC 

The point where the first order beam corresponding to the 
minimum passband frequency intersects the zero order dif- 
Fraction pattern determines a minimum detectable signal level 
due to zero order interference effects. This minimum signal 
value relative to the axis intensity of the zero order beam 
was calculated as follows for the NPS a-o spectrum analyzer 


assuming an octave bandwidth: 


ER 
I =I (98_._) = minimum detectable signal 
S o min 
Then 
I, 
T (o) -44 dB 
O 
where: 
W= D9 ._/A -= 51 
min’ o 
D = 7 mm 
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See. = A f / (2V) 


min Omin 
A u en 
f = 80 MHz 

O 


V = 3.63 km/sec 


It is important to note that the above calculated figure is 
independent of wavelength. The only way to lower minimum 
detectable signal level due to diffraction is by sidelobe 
reduction. For systems in which sidelobes are sufficiently 
attenuated, the minimum required signal would be limited by 
Brillouin scattering. The scattered power due to this effect 


is given in the following equation / 28 /: 


O 
ENEE 
Pr = total power scattered in all directions 
E = incident power 
F = scattering EE proportional to T? 
T = temperature in degrees Kelvin 
L = interaction medium length 


In the subsequent calculation of minimum ratio, F was not 
available for Lead Molybdate so a representative coefficient 


for a very strong scatterer was used. 


tg 


_ 2FLS 


= Es = normalized part of incident power delivered 


O 


into one resolution cell (54) 
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~ -130 dB 


where: 
F= 40 om? 
L= ./ cm 
S = Ab(1 - cosA9) 
AQ = A¿/E 
ro = „6323 um 


If the two power and intensity ratios are taken to be essen- 
tially equivalent measures, then the zero order diffraction 
sidelobe level was clearly much stronger than the noise level 
due to Brillouin scattering. The upper linearity limit was 
determined by the maximum single frequency input allowable 

in a multifrequency environment which has been shown by Hecht 


/ 24 / to be: 


LAS al — 20 dB 


The resulting dynamic range (DR) calculated for the NPS Lab 
deflector with this condition was then: 


DR = 24 dB 


In general, it would appear from the above typical calcu- 
lations that systems, which do not employ special sidelobe 
attenuation techniques, would be zero order diffraction 
limited in dynamic range. However, F is very temperature de- 
pendent so a unit might begin operation limited by zero order 
sidelobes and gradually suffer a reduction in dynamic range 


due to device heating. 





With discussion of diffraction efficiency now completed, 
concentration is shifted to information capacity. This per- 
formance measure is defined as the number of beam spots 


resolvable in the image plane and is determined by: 


20 = A/A = deflection angle from zero order beam 
SE" di p in 
= À 8 E/V 
N = A¢ọ/49 
= Aft = number of resolvable beam spots in image plane 
where: 
t = D/V 


D = effective aperture width 

a= acoustic velocity 
eer ge information capacity is desirable but difficult to 
Seran, as only two factors affect its value. Bandwidth 
improvement considerations have already been discussed at 
length under diffraction efficiency and so are not covered 
Mere. The length of time t required for a wavefront to travel 
the effective length of the sound column is the principal 
HEEL of this section. To lengthen tT, a low acoustic velocity 
would apparently assist matters, but acoustooptic materials 
which have low sound velocities also tend to have high acoustic 
Mess factors. This attenuation limits the deflector's effec- 
tive aperture width and gives the sound column an exponential 


amplitude profile shown by: 


A(x) = a (55) 





where: 


A(x) = normalized acoustic amplitude 
X = distance from transducer along sound column 
oO = acoustic loss factor 


= r£°v/{20109, ge) nepers/sec 
T = material acoustic loss constant (dB/{em-GHz~}) 


V = acoustic velocity 


In order to gain a rough measure of the effective aperture's 
relation to acoustic loss factor, an arbitrary width was 
chosen where A(D) equals 2/e or approximately 1//2. Brief 


algebraic manipulation yields: 


v, 
N 


{Vlog_2}/a (56) 


Pt 


The number of resolvable frequencies may now be rewritten as: 


> 
N 


are 
O O 


ona) 


K/S 


where: 


bx 
il 


llL_/{VILÌ 


This result clearly shows better acoustic materials have 
smaller VI products. It also points out for a given medium, 
maximum information capacity is inversely proportional to 
center frequency. The implication for operating systems is 
that once the "best" material has been selected for a normal 


diffraction Bragg cell then only a reciprocal relationship 
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exists between the system's surveillance bandwidth and infor- 
mation capacity. One must be sacrificed by the same percentage 
the other is improved. However, this last conclusion repre- 
sents the ideal case, inferior performance would be observed 
for systems which allow thermal gradients to form during 
periods of high power input. In Ref. 10, Uchida presented 
findings on this subject which showed in one case as much as 
Aye percent drop in information capacity for an input power 
increase of six percent. Proper heat sinking was reported 
to significantly reduce this heat dependence. 

Quick response is the last general topic in the discussion 
of Bragg cell performance factors. This performance measure 
is desirable in some applications, such as Bragg cells intended 
to be beam modulators rather than deflectors. On the other 
hand, high resolution surveillance spectrum analyzers would 
not be noted for response agility as information capacity would 
directly suffer. A tradeoff must be made for optical excisors 
as both resolution and response sveed are desirable here. For 
END ur pose Of this work, the sharpness of Bragg cell response 
is defined as one over T --the rate at which single resolvable 
frequency inputs may be "hopped" and still attain maximum in- 
Bau vyadirtraetion ratio. It is a measure of how 'instan- 
taneous" the displayed spectrum is. The diffraction intensity 


of short pulses is determined by / 29 /: 


D 
TEO) - SARE) I (x) dx (58) 
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where: 


I_(t) = "Short Pulse" instantaneous light intensity 
response 
A(x,t) = Q lative diffraction Of signal In medium as a 


EE ion of transited distance and Cime 


1, (x) incident light beam cross section in aperture 


D 


effective aperture width 


Short pulse attenuation effects were particularly noticeable 
when observing the spectrum of a sinusoidally modulated FM 
imput signal shown in Figure 14. The input power level re- 
quired an increase of nearly two orders of magnitude in 
contrast to the power needed for a single beam. For a quali- 
tative overall performance comparison, a display of the same 
Signal spectrum produced by a conventional scanning super- 
heterodyne spectrum analyzer is presented in Figure 15. 

In this section, a thorough compilation and examination 
of all important parameters was attempted. To summarize the 
findings: large information bandwidth, quick response, high 
diffraction efficiency and large dynamic range are functions 


of geometry, material constants and tradeoffs. 


Dewey PHOTODETECTOR PERFORMANCE 

The Fairchild @CDI21H and Reticon CEPD=1/28 photodetector 
arrays were the principal detectors examined in this work for 
spectrum analysis application. The two 1728 photosite arrays 
have many similar features of interest, especially in their 
internal working structures. Each detector photosite of both 


arrays integrates the current generated by the incident light. 
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Figure 14. NPS A-O Spectrum Analyzer 
Display of Wideband FM Signal 





Figure 15. Conventional Spectrum Analyzer 
Display of Same Wideband FM 
Eeer 
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For this charge storage, the Reticon array photosites use 
photodiode junction capacitance while the Fairchild device 
elements rely on dielectric capacitance. Both arrays transfer 
the stored charge to CCD analog shift registers for further 
transfer to amplifying circuits and eventual display. This is 
accomplished by dividing the array elements into even and odd 
groupings. Even elements transfer to a shift register on one 
Side of the photosite array and odd ones transfer to an 
Opposite register. The first difference in operation occurs 
When the registers are transferring data off chip. The CCD 
chip interleaves the even and odd charge packets prior to off 
chip exit. The CCPD ship does not, but rather leaves this 
"EE cOr external circuitry. Diagrams of both chip operations 
are shown in Figures 16 and 17. 

Examination of respective data sheets for the two devices 
reveals near equivalence in all specifications save two. The 
CCD array was a denser and faster chip. It required fourteen 
percent less area per photosite and used rectangular vice 
square shaped elements with the short dimension along the 
array face. This translated into a twenty percent resolution 
improvement over the photodiode chip. Also, the maximum output 
data rate of the Fairchild device was twice that of the Reticon 
chip. This meant better resolution of input signal temporal 
characteristics. 

Dynamic range for both chips was defined as the individual 
array element's ratio of saturated output voltage to background 
noise equivalent voltage. The advertised dynamic range for 


Sach device was given as 27 dB. This figure was verified by 
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Figure 16. Fairchild Charge Coupled Device 
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Figure 17. CCPD Photodetector Array Chip Schematic 
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two separate methods for the NPS spectrum analyzer CCD 
detector. The first method utilized the calibration marks 

on the rf input signal generator for this measurement. The 
subsequently observed dynamic range was 27.5 dB plus or minus 
one dB. The second technique employed neutral density 
filters to attenuate the saturated output signal down to the 
noise level and yielded a 26 dB range. However, slight beam 
spreading due to filter presence was noted which further 
attenuated the signal and accounted for the lower assessment. 

Brief experience was gained with a smaller Reticon CCPD- 
Peewecdip and related support circuitry before working with 
the Fairchild CCD-1728 chip and development board. The photo- 
diode array required slightly less effort to fine tune than 
the Fairchild model, however this may have been due to the 
simoler nature of the smaller device package. The larger 
Fairchild unit was used extensively over an eight month period 
and performed with 100 percent reliability throughout this 
time. The author believes either device technology would 
perform commendably in spectrum analyzer applications, but is 
biased in favor of the straight CCD approach in view of the 
above noted advantages. 

RE Ee neted in he introduetion s list of 
specific goals, require a detector much different from the 
devices discussed above. Most excisor applications require 
wide information bandwidths of the kind offered by photomul- 
tiplier tubes (PMT) available on an "off the shelf" basis. 

For example, Varian markets typical PMTs which have bandwidths 


of a GigaHertz and dynamic ranges of 114 dB. These specifications 
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mean that a given system's performance will be most likely 
limited by the Bragg cell dynamic range and bandwidth. With 
these comments, focus is now placed on performance modeling 


Ee optical excision process. 


E. ELECTRONIC AND OPTICAL EXCISOR PERFORMANCE MODELS 

Optical excision performance is measured in terms of two 
gauges which are: how well offending signal frequency com- 
Pemencts are removed from the output, and how much the original 
signal of interest is corrupted by the loss of frequency com- 
ponents shared with the offending signal. Both gauges are 
Fesmaraly functions of the electronic or optical clipper 
selected for use in the particular system. In this analysis, 
electronic clippers are assumed to completely remove inter- 
ference signals chosen to be excised since nearly infinite 
attenuation would be available at the output ports. This 
same claim cannot be made for optical clippers because of 
diffraction limitations. The narrow obstacle of the optical 
clipper impedes the progress of the intended beam in much the 
same manner as a tree or telephone pole positioned between 
two people impedes dialogue. The observed diffraction effects 
are advantageous for conversation but not for excision. 

Examination of optical clipper efficiency was approached 
by assuming the obstruction was much narrower than the im- 
io ino diffracted beam. In the limit, this situation may 
be described as a narrow obstacle in a coherent plane wave 
light field. Clipper effectiveness then becomes a function 


of the narrow obstacle width and the distance to the image 
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plane depicted in Figure 18. The resulting image plane in- 
tensity diffraction pattern provides a visual measure of 
en pper utility. This diffraction pattern was calculation in 
terms of Fresnel integrals where the narrow obstacle was 
considered to be the intersection of two overlapping infinite 
half planes. For further information, brief developmental 
guidance to derive the following image plane intensity profile 
is given in Ref. 17. The normalized light intensity profile, 
I(U) , along the u axis in the image plane (Fig. 19) is given 
by: 

Me nad any A (59) 
where: 


u 
a, li I ea dw 
=2 


u 
B= J d re? dw 
EZ 


Z as 
oo 73 i) 
O 


= normalized integral limit 


y >= actual integral limit mapped to U, 
Za y~ actual obstacle width 
r = distance from obstacle to image plane 


A = light wavelength 


This equation was approximated by Simpson's rule and used 
in a Fortran program to plot image plane intensity profiles 
shown in Figure 20 through 26. The figure sequence shows the 
effect of an obstacle which is moved steadily away from the 


image plane until its effect is negligible. The program 
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Figure 18. Narrow Obstacle in Plane Wave Light Field 
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Figure 19. Narrow Obstacle Diffraction Model Schematic 
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used for plotting is listed in Appendix C. 

To analyze the effects of signal corruption by the excisor, 
the reader is directed to consult Ref. 30. In this work, 
Sussman and Ferrari conducted a study of Direct Sequence 
Spread Spectrum correlation peak attenuation due to the pre- 
sence of an inband notch filter. They found that interference 
frequency spikes, which have bandwidths of up to three percent 
of the spread spectrum signal bandwidth, can be excised with 
only one dB drop in the output correlation peak. This worst 
case result was obtained for a notch filter centered on the 
spread spectrum signal center frequency. Much proportionately 
wider excision bandwidths were tolerated for slightly offset 
notch filter center frequencies. In conclusion, their analysis 
indicated the comparatively slight signal degradation due to 
the notch filter was much preferable to the degradation caused 


by the narrowband interference. 


ER 








III. SUMMARY/RECOMMENDATIONS 


A. RESULT SUMMARY 

Laser beam propagation, Bragg cell performance, photo- 
detector performance, and optical excision models were the 
main topics addressed in this work. For the laser, optimal 
light beam profile behavior in the image plane was found to 
be a function of both Gaussian beam propagation and aperture 
truncation effects. Observations of an assembled spectrum 
analyzer were made which compared closely with calculated 
values. And, two methods of beam sidelobe attenuation were 
proposed for signal resolution improvement. Bragg cell per- 
formance was measured in terms of diffraction efficiency, 
information capacity, and response agility. Also, device 
dynamic range was calculated for a multifrequency input signal 
environment. Additionally, beam steering and diffraction 
efficiency effects of light-sound beam geometric factors were 
studied. Photodetector array performance of charge-coupled and 
photodiode technologies was compared for spectrum analyzer 
application. Optical excisor models were discussed in terms 
of interference removal effectiveness and subsequent corruption 
of residual information bearing signals. A Fresnel approximation 
model of optical clipping was developed and programmed for com- 
puter plot generation. And finally, signal degradation was dis- 
cussed for DSSS signal correlation peak attenuation due to 


excision. 
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B. RECOMMENDATIONS 

The following remarks are specific recommendations for 
optimal performance of spectrum analysis, and optical excision 
of narrowband interference. To reduce zero order sidelobe 
interference, and thereby increase dynamic range, use of Bragg 
cell apertures which do not mabruptly change in optical trans- 
mittance at the window edges is recommended. Choice of a 
deflector interaction medium having high refractive index, 
large elastooptic coefficient, small mass density, and small 
acoustic velocity-material acoustic loss constant produce is 
also desirable. Either photodetector technology was found 
acceptable with a slight advantage held by the charge-coupled 
device array. Unfortunately, no guidance can be given with 
regard to the choice of a "best" clipper for optical excision 


as developments in this area are still progressing quickly. 


Ge AREAS FOR FURTHER» ANALYSIS 

Further effort related to this thesis is required in 
Gaussian beam truncation analysis and electronic/optical 
clipping model development. A Discrete Fourier Transform 
(DFT) computer model of aperture windowing of Gaussian beams 
1 ne possible way of promise to study truncation effects as 
Hat ic solutions are exceptionally difficult to arrive at. 
Determination of window functions which "best" attenuate 
sidelobes and maximize information capacity would be highly 
desirable. Much further work also remains in optical excision 
modeling. Another DFT computer simulation could be used to 


render a Fraunhofer apploximation model of obstacle-light beam 
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MI. aCtion in order to further clarify optical clipper 
effectiveness. This approach could easily handle the case 

where the obstacle width would be equal to the incident light 
beam width. Lastly, a similar model could be used to determine 
nearly exact degradation effects of the excision process. These 
results could then determine optimal utilization of the exci- 


sion process. 


De CONCLUSION 

The performance analysis guidelines developed in the body 
of this thesis provide a realistic view of acoustooptic poten- 
tial for spectrum analysis and interference excision applica- 
tions. Also, several acoustooptic spectrum analyzers have 
been made with commercially available components and readily 
provide wideband, high resolution coverage not easily obtained 
in Conventional systems. Lastly, work in optical excision 
techniques has progressed substantially and ongoing efforts 
continue to show performance improvement. In view of these 
facts, further research effort in this area is very likely to 
provide many more valuable benefits for applications involving 


wideband, high density signal environments. 


OR 








APPENDIX A 


BEAM PHASE FRONT RADIUS AND SPOT SIZE PLOT PROGRAM 


This Fortran program plots Gaussian beam phase front 
radius and spot size as a function of distance from beam 


waist. Selectable parameters are: 


LMBDA = laser wavelength 


NPTS = number of plot points desired 


WZERO waist beam spot size 


DCRM distance increment 


Prospective users at the Naval Postgraduate School should 
refer to the W. R. Church Computer Center's Versatec plotter 
routine user manual for required job control cards and plot 


dimensions. 


LO 
LO 








100 


399 


BEAM RADIUS AND WAIST AS A FUNCTION OF DISTANCE 


DIMENSION R(500), X(500), W(500) 
INTEGER*4 ITB(12)/12*0/ 
REAL*4  RTB(28) /28*0.0/ 


RPTS = 500 

REAL LMBDA 

LMBDA = .6328 E-06 
WZERO = „44 E-03 


EE + I 
DCRM = 0.03 
PI = 3.141592654 


Te 

Z(1) = 1.0 E-03 

R(I) = Z(I) + (PI * WZERO**2/LMBDA) **2/Z (I) 

A = SOTR(1.0 + (LMBDA * Z(I)/(PI*WZERO**2))**2) 
W(I) = WZERO * 2 

ee 


TEO TT EOC T EMT) CGO TO 999 
Z(I) = 2(1 - 1) + DCRM 


ITB(3) = 8 
ITB(4) = 4 
RTB(1) = "5 
RTB 2 ome 5 E-03 
GO TO 100 


CALL DRAWP (NPTS,Z,W,ITB,RTB) 
EB (232-5570 

CALL DRAWP (NPTS,Z,R,ITB,RTB) 
STOP 

END 


100 








APPENDIX B 
GAUSSIAN PROPAGATION PROGRAM 


This TI-58, TI-59 hand-calculator program calculates 
Gaussian beam phase front radius and spot size for propaga- 
tion through any number of thin lenses spaced at arbitrary 
distances. Selectable parameters are stored in the following 


registers: 


*R00 = laser wavelength 
ROl = beam spot size 
ROS = lens focal length 


RO4 = distance increment 


Intermediate results are stored in: 


R06 real part of complex q 


R07 = imaginary part of complex q 
* Common units are used throughout the program 


OPERATING SEQUENCE 

The wavelength of the laser is the first necessary entry. 
This value is entered into program memory by loading it into 
the display register and pressing label A. The initial beam 
spot size and phase front radius are entered by first keying 
spot size into the display and pressing B. Immediately after 
this, the phase front radius is loaded in the display and 
entered by pressing the R/S key. If the user desired to start 


at the waist then an infinite phase front radius can be 


10% 








entered by loading any number greater than or equal to 1000, 
as all such values for phase front radius are treated as 
Wisen te, 

The above portion of the sequence determines the initial 
conditions for beam propagation. From this starting point, 
any amount of propagation distance may be entered by loading 
the additional path length desired in the display and pressing 
D. The newly calculated beam spot size is stored in the 
display register and the corresponding phase front radius is 
stored in the t-register. Propagation through thin lenses is 
handled in a similar manner by entering the focal length in 
the display and pressing C. Spot size and phase front radius 
are read out in the same manner as for additional propagation 
distance. 

The above comments complete the operating sequence des- 
Fi tion. This is now followed by a listing of the program 


steps by their functional sections. 


SEQUENCE FOR ENTERING LASER WAVELENGTH: 
*Lbl A STO 00 R/S 
SEQUENCE FOR ENTERING INITIAL BEAM SPOT SIZE AND PHASE FRONT 
RADIUS: 


le STO 01 R/S soe 02 T000 x<- E RCL 02 *x>t INV 


Men O02 1/x + RCL 02 = Mid enema we RED Ol 
x°)) x*)) STO 85 

Ed = *m + RCL OI yl + (ReBEO2 MI < SE (RCL 00 
+ *T + RCL Ol a Ta STO 06 


w02 








een) 4<-t CLR R/S 


*Lbl INV CLR STO 05 (RCL 01 x* o nz RCL 00) 


STO 06 CLR X<>t CLR R/S 


SEQUENCE FOR PROPAGATION THROUGH SPATIAL DISTANCE: 


DLD + RCL 05 = STO 05 E 


SEOUENCE FOR CALCULATING BEAM SPOT SIZE AND PHASE FRONT 
RADIUS FROM COMPLEX "q": 


*Lbl E ((RCL 05 xÉ + RCL 06 x°) + RCL 05) STO 02 


ME: (RCL 00 + *m + RCL 06 x (RCL 05 x° + RCL 06 =a 


ee" sTo 01 R/S 


SEQUENCE FOR TRANSMISSION THROUGH THIN LENS: 


MoL C STO 03- (RCL 03 wi x RCL 06 + ((RCL03 - 
ECE 05) wi + RCL 06 x? ire ORC xa) REE OS 
Meer « RCL 05 x (RCL 03 = RCL 05) = RCL 03 
A RCL 06 SCH + (REL OIE BOB O05) = + RCL 06 =. 


SED OS x<>t STO 06 E 
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APPENDIA C 


BANDSHAPE PROGRAM 


This Fortran program plots Bragg cell passband shapes for 
arbitrary integer ratios of interaction length (L) to charac- 
teristic interaction length (Lo). Selectable parameters are: 

NPTS = number of plot points desired 


SDLTAX normalized frequency increment 


id = normal'ized start frequency for plot 


The number of plot lines is controlled by the limit of the 
second DO loop. Prospective users at the Naval Postgraduate 
pool should referto the W. R. Chruch Computer Center's 
Versatec plotter routine user manual for required job control 


cards and plot dimensions. 


* normalized to 1 
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ee SINXTSOR OVER X SOR 


10 


100 


200 


DIMENSION X(0201), F(0201) 
INTEGER*4 ITB(12)/12*0/ 

REAL*4 RTB(28)/28*0.0/ 
EQUIVALENCE (TITLE,RTB(5)) 
REAL*8 TITLE(12)/'BAND SHAPE'/ 


NPTX = 201 

DLTAX = .0L1 

PI = 30141592654 
lee 

QPRME = 2.0 PL 
Q= 0.0 


Dom I Mol 201 

X(I) =X(I -1) + DLTAX 
CONTINUE 
DO 200 J 
Q = QPRME 
Donner. = 1,2203 
E ER 
ee Se), 

Hat) EO C/B Vase 
CONTINUE 

E EE ESRI 
EE EES 
TEO O EE 
BABS) 
ITB (4) 
RER 
ITB(10) = 
Fre) 
EENS 
CALL DRAWP (NPTS,X,F,ITB,RTB) 
CONTINUE 

STOP 

END 


, 6 


+ Il 
OH 


zelt 
WNW 


8 
i 
2 


2 
.4 
.4 


105 





APPENDIX D 


NARROW OBSTACLE DIFFRACTION PROGRAM 


This Fortran program plots diffraction patterns caused 


by a narrow 


obstacle in a plane-wave light field using a 


Simpson's rule approximation algorithm. Selectable parameters 


are: 


LMBDA = 


YMIN 


DLTAY 


RZERO = 


INPP = 


NEP 


Prospective 


laser wavelength 

left plottimg boundary limit 

obstacle width 

distance from obstacle to image plane 
total number of plotting points 


total number of integration points 


users at the Naval Postgraduate School should 


refer to the W. R. Church Computer Center's Versatec plotter 


routine user manual for required job control cards and plot 


dimensions. 
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C NARROW OBSTACLE DIFFRACTION PROGRAM 


200 


100 


300 


DIMENSION Y(500), X(500), X(500) 

INTEGER*4 ITB(12)/12*0/ 

REAL*4  RTB(28)/28*0.0/ 

EQUIVALENCE (TITLE,RTB(5)) 

REAL*8 TITLE(12)/'DLTAY EQ 2 MICRONS RZERO EQ 19 MM'/ 
REAL LMBDA 


LMBDA = .6328 E-06 
YMIN = -40.0 E-06 
DLTAY = 2.0 E-05 


RZERO = 19.0 E-03 

HPI = 1.570796327 

TNPP 500.0 

TNIP 8000.0 

STR = SORT(2.0 / (LMBDA*RZERO) 
UMIN = YMIN * STR 

DEAU SIDILTAY * STR 

DIT DLTAU / TNIP 

DPP = -2.0 * UMIN / TNPP 


Tel 

TUMIN = UMIN 
PUMIN = UMIN 
SUMI = 0.0 
SUM2 = 0.0 


CITS TCOS (HPI (TUMIN** 2) ) 

C2 = GOs (oP l= (LUMIN+DII) **2) 

EE Ee EE ee) 

SUMIL = el + 4.0*C2 + C3) *DII/3 + SUML 

EEN + 2.0*DII 

EN EIE EE. (PUMIN + DLTAU [= 2.0*DII)) GO TO 100 
TOMINI RUMIN 


SUS SIN API" (TUMIN**2)) 
SMS UNTAR TS (TUMIN+DII)**2) 
SON SAP Es (TUMIN+2.0*D1TI)**2) 


SUM2 = (Sl + 4.0*S2 +S3)*DII/3 + SUM2 

TUMIN = TUMIN + 2.0%DII 

IF (TUMIN LE, (PUMIN + DLTAU - 2.0*DII)) GO TO 300 
dT EG 251) GO TO 999 

X(I) = (SUM2-SUMl)**2 + (2.0 - SUM1 -SUM2) **2 

Z(I) = X(I) / 4.0 

ASSI = Z(I) 

AR PUMIN + DLTAY / 2.0) / STR 


TONE) — YT) 
I=I+ 1 

PUMIN = PUMIN + DPP 
NPTS = 500 

ITB(3) = 8 

ITB(4) = 4 

A 22 

ITB(10) 


= 2 
ITB(11) = 4 
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RTB(1) = 1.0 E-05 


RTB(2) = .5 
GO TO 200 

999 CALL DRAWP (NPTS,Y,Z,ITB,RTB) 
STOP 


END 





Or, 


TIE, 


Im). 
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